Spindles and active vortices in a model of confined filament-motor mixtures 
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Robust self-organization of subcellular structures is a key principle governing the dynamics and 
evolution of cellular life. In fission yeast cells undergoing division, the mitotic spindle spontaneously 
emerges from the interaction of microtubules, motor proteins and the confining cell walls, and 
asters and vortices have been observed to self-assemble in quasi-two dimensional microtubule-kinesin 
assays. Their is no clear microscopic picture of the role of the active motors driving this pattern 
formation, and the relevance of continuum modeling to filament-scale structures remains uncertain. 
Here we present results of numerical simulations of a discrete filament-motor protein model confined 
to a pressurised cylindrical box. Stable spindles, nematic configurations, asters and high-density 
semi-asters spontaneously emerge, the latter pair having also been observed in cytosol confined 
within emulsion droplets. State diagrams are presented delineating each stationary state as the 
pressure, motor speed and motor density are varied. We further highlight a parameter regime 
where vortices form exhibiting collective rotation of all filaments, but have a finite life-time before 
contracting to a semi-aster. Quantifying the distribution of life-times suggests this contraction is a 
Poisson process. Equivalent systems with fixed volume exhibit persistent vortices with stochastic 
switching in the direction of rotation, with switching times obeying similar statistics to contraction 
times in pressurised systems. Furthermore, we show that increasing the detachment rate of motors 
from filament plus-ends can both destroy vortices and turn some asters into vortices. Based on our 
findings we argue the need for a deeper understanding of the microscopic activities underpinning 
macroscopic self-organization in active gels and urge further experiments to help bridge these lengths. 
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I. INTRODUCTION 



Filamentous proteins are prevalent within eukaryotic 
cells and perform a variety of crucial tasks relating to cel- 
lular integrity, locomotion, transport and division [HE]- 
Such tasks are often active in that they can only proceed 
in concert with energy-consuming mechanisms, including 
directed filament growth and motor protein-generated 
tension, placing such processes outside the realm of equi- 
librium thermodynamics [3j. Self-organisation of motor 
protein-filament mixtures will be selected for when it ro- 
bustly reproduces static or dynamic structures beneficial 
to the cell's viability. An example is the mitotic spin- 
dle that forms during division of fission yeast cells. It 
has been shown that this bipolar structure, consisting 
of microtubules emanating from spindle pole bodies to- 
wards an overlapping midplane region, exists and func- 
tions essentially as normal even in cells with no nucleus- 
associated microtubule organizing center [HIS]. The plau- 
sible conclusion is that the interaction between filaments 
and motor proteins in the confined cell geometry controls 
the location of the pole bodies. For budding yeast this 
self-organisation scenario has been reinforced by the evo- 
lution of more sophisticated regulatory mechanisms [6], 
and egg cell extracts from the amphibious genus Xeno- 
pus can generate a well-formed spindle apparatus despite 
entirely lacking cell walls |7 9 . Nonetheless an under- 
standing of the principles underlying self-organization of 
bioflaments driven by motor proteins in confined spaces 



is of direct relevance to many organisms |10| . 

Given the complexity of real cells it is often advan- 
tageous to consider simplified model systems, and this 
approach has been adopted to investigate the role of 
confinement in filament-motor mixtures. Experiments 
on growing microtubules confined to spherical emulsion 
droplets revealed a droplet-size dependency on the ob- 
served structure [TT]: Droplets larger than f« 29/xm in 
diameter contained asters with the polar microtubules 
pointing towards the centre, controlled by the motor pro- 
tein dynein, whereas smaller droplets were found to con- 
tain semi-asters with the aster's focus near the interface. 
These findings demonstrate that the degree of confine- 
ment can partly determine structure formation, but as 
motor density and speed were not control variables in 
these experiments their infiuence could not be assayed. 

A strikingly non-equilibrium property of filament- 
motor mixtures is their ability to spontaneously gen- 
erate flows due to their active components, even in 
the absence of boundary driving forces [I2l414j . Assays 
of microtubule-oligomeric kinesin mixtures in a quasi- 
two dimensional geometry with fiat, parallel confining 
walls found a dynamic rotating structure denoted a vor- 
tex \15\ I16j . Accompanying simulations of semiflexi- 
ble filaments |16| and subsequent hydrodynamic theo- 
ries [T71 [TB] appeared to reproduce the observed struc- 
tures. However, as discussed in Rcf. [19], it is unlikely 
that the simulations of Surrey et al. [16] and the theo- 
ries and simulations of Ref. [TT] [TH] describe the same 



type of vortex, because the hydrodynamic theories are 
based on a nematic order-parameter description, while 
simulations of semi-flexible filaments in Ref. [TB] neglect 
self-avoidance (and thus nematic order). Simulations of 
self-avoiding filaments strictly in two dimensions showed 
no evidence of a vortex state [19 . The microscopic 
picture underlying vortex formation thus remains un- 
known. Gliding assays of filaments along motor beds 
permit quantitative comparison to models HD] and at 
high concentration exhibit vortex-like 'swirls' (SUES], al- 
though in this situation the active forces are unbalanced 
monopoles, unlike dipoles generated by motors connect- 
ing two filaments in the bulk [53j. Vortex-like motion 
is often observed in self-propelled systems such as bac- 
terial swimmers |24H26) . but with differing microscopic 
mechanisms. 

It is apparent that the combined influence of con- 
finement and activity on structure formation and spon- 
taneous fiows in filament-motor mixtures is presently 
not well understood. Our aim here is to acquire a 
deeper understanding of this problem in a broad sense, 
not restricted to any one biological realisation, i.e. 
microtubule-dynein or actin-myosin. It is therefore de- 
sirable to study model systems in which all parameters 
can be freely varied. The application of continuum equa- 
tions, which are coarse-grained over lengths much larger 
than the filament length L, to structures of only a few 
L in spatial extent is not guaranteed to be successful. 
We therefore adopt a discrete numerical model in which 
motors and filament segments are explicitly represented, 
and all physical mechanisms that are potentially rele- 
vant (steric hinderance, thermal fluctuations etc.) are 
incorporated. This model is an extension of one pre- 
viously employed in two dimensions |19j . where it was 
found to produce some signatures of active gels such as 
super-diffusion and anomalous small wavelength density 
fluctuations, but not vortices. 

We consider arrays of filaments confined to a quasi- 
two dimensional cylinder, with a height of a few fila- 
ment diameters which permits filament overlap, and an 
external pressure at the curved walls. We then system- 
atically vary the motor density, speed and applied pres- 
sure. Four steady-state configurations arise within the 
covered parameter space, including an aster and semi- 
aster as observed in confined emulsion droplets [llj . and 
also a spindle-like state that spontaneously emerges from 
the motor-filament interaction in the confined geometry, 
possibly reproducing the fission yeast observations [H [5] . 
These states are described in Sec. |III A| along with a 
fourth nematic state that links to known equilibrium 
phases. We also find a fifth, vortex state associated with 
a definite rotation of filaments about a fixed center that 
appears to be always transient. The existence and prop- 
erties of these vortices are characterised in Sec. IIIIBI To 
highlight the important role played by motors at filament 
plus-ends, we independently vary the detachment rate of 
motors from plus-ends in Sec. |IV| and show that vorticity 
is associated with a critical fraction of plus-ended mo- 



tors. The observation of vortices in fixed volume systems 
described in Sec. Iv] confirm that they are driven at least 
partly by motor motion and not boundary fluctuations, 
and in Sec.|VI|we discuss possible future directions. 



II. MODEL 

We consider a system of N semiflexible polar fllaments, 
which can be connected by motor proteins. Each fila- 
ment consists of M = 30 monomers separated by a bond 
length b with Hookean bond potentials with a spring 
constant 100 k-gT /h. Self-avoidance of filaments is in- 
troduced by repulsive Lennard- Jones potentials with di- 
ameter a and energy parameter e = 5 k^T. A natural 
choice is (J = 6. Semi-flexibility is described by curvature 
elasticity with bending rigidity k, = 200bkBT such that 
the persistence length £p = n/ksT is ip = 20L/3 with 
L — Mb the filament length. Motors are modeled as two- 
headed Hookean springs with a spring constant k^T/b^ 
and dynamics defined by four rates as shown in Fig.nja): 
(i) The attachment rate /ca for a motor to attach to two 
monomers within a predefined range; (ii) the detachment 
rate ^d of each head independently from its filament (de- 
tachment of either head results in removal of the whole 
motor from the system); (iii) the movement rate fcM of 
each head independently towards the filament's [-|-]-end, 
and (iv) the detachment rate /ce for motor heads already 
at a [+] end. The movement rate is attenuated by an ex- 
ponential factor e-^^/f^BT .^^j^jj /s^jt; h^q change in motor 
spring energy for the trial move. Except where otherwise 
stated, /ce = ^d below. 

The quasi-2D simulation cell has parallel confining 
walls normal to the z-axis spaced 56 <C i apart, restrict- 
ing filament orientation to be approximately in the x~y 
plane while permitting a degree of overlap; see Figs, mb) 
and (c). These walls repel the monomers with the same 
Lennard- Jones non-bonding potential as for filaments. 
N — 175 filaments are placed in a radial aster config- 
uration with all [+]-cnds pointing towards the center, 
in three parallel layers with roughly 66-68 filaments per 
layer (note there is some stochasticity in the initial con- 
ditions). This initial condition was chosen to promote 
the formation of asters and vortices, but does not in- 
hibit other structures as described below. This initial 
array is surrounded by an elastic wall with both bend- 
ing and stretching energies, with the elastic coefficients 
chosen to ensure a circular shape is approximately main- 
tained throughout. The stiffness with respect to changes 
in wall area is chosen to be sufficiently small such that, 
for the parameters considered below, wall retraction is 
countered dominantly by filaments repulsion and not the 
intrinsic wall elasticity. It does however remove the pos- 
sibility of wall nodes merging and causing a numerical 
divergence. Stretching and bending energies were de- 
fined similar to above with monomers replaced by wall 
nodes, with stretching and bending coefficients k-^T/l'^ 



and 2000 A:bT resp., with 



nb the initial separation 
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FIG. 1: (a) Summary of key model parameters including the 
rates of motor attachment /ca and detachment ku, and the 
bare stepping rate ku- See text for details, (b) Plan view 
showing the filaments oriented with their light-shaded [+]- 
ends towards the center. The arrows denote the external 
pressure acting on the circular elastic wall, (c) Side view 
of the same, showing the confining walls perpendicular to the 
z-axis. 



between nodes. 

The filaments, motors and elastic wall are all updated 
stochastically. The filaments obey Brownian dynam- 
ics I27j governed by an effective monomer friction co- 
efficient 7 for hydrodynamically anisotropic slender ele- 
ments, i.e. with a 2:1 ratio between implicit solvent drag 
perpendicular to the filament axis (= 27) to the par- 
allel direction (= 7) [28]. Wall nodes move by Monte 
Carlo Metropolis moves applied to the {x, y) coordi- 
nates of 80 nodes initially equispaced along its contour. 
The energy for these moves includes the wall elastic and 
wall-filament interaction energies, and a pressure- volume 
term PV where P > is a fixed parameter for each 
run. To check for convergence with time, various scalar 
quantities, such as the number of motors per filament, 
were checked to be constant within noise when plotted 
against log(t). In addition, the mean squared rotation 
{Mf{t,t + At) = iV-i E^i Mt + At) - 0,{t)\^, with 
di the angle between filament Vs centre-of-mass and the 
nominal centre of the box (more precisely, the mean of 
all wall nodes) relative to some fixed axis, was checked to 
no longer to vary with t to within noise. Stationarity was 
not achieved for the vortex states, for which alternative 



III. CONTROLLED PRESSURE WITH Ae = ku 

Results are presented here in terms of the normalised 
attachment rate kx/kjy, the normalised motor rate kyiTb 
where r;, = LJry/AkBT is the approximate time for a fil- 
ament to freely diffuse one monomer distance; the nor- 
malised pressure P/Pq with Pq = e/b^ with e = SkeT the 
Lennard- Jones repulsion energy, and, where relevant, the 
scaled end-detach rate kE/ku. We first describe the ob- 
served steady-state configurations before turning to con- 
sider the transient vortices. 



A. Stationary states 

For the parameter space considered, we observe four 
classes of steady-state configuration as shown in Fig. [2] 
For a low density of fast motors, spindles are observed as 
in Fig.^a), which crossover to a radially-symmetric aster 
as the motor density is increased as show in Fig.l2jb). For 
slower motors, we observe a nematic at low motor den- 
sities and semi-asters at high motor densities, as shown 
in Figs. l2jc) and (d) resp. Semi-asters typically arise for 
higher pressures than asters and are more compressed, 
consistent with the emulsion experiments of Pinot et 
al. [11 and justifying our use of the term. Movies demon- 
strating the spontaneous emergence of all of these states 
from the initial conditions are provided in the supplemen- 
tary information |29| , for exactly the same parameters as 
in Fig. [I 

To quantify to which state a system belongs, each fil- 
ament's polarity vector is projected onto the x-y plane 
to give a two-dimensional unit vector aligned towards the 
[-t-]-end. This is averaged over all filaments whose centers 
of mass have azimuthal angle 6 with respect to the center 
of the system, giving rise to the mean orientation p{9). 
This is then decomposed into angular mode vectors a^ 
and b„i. 
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from which can be defined the mode amplitudes Q„ 



Qm ~ n 2 ^^rn 



+ 6™)- 



(1) 



(2) 



measures were employed as described in Sec. IIIB 



The Qm are invariant under global rotations of the whole 
box. 

To determine the corresponding state, the measured 
Qni up to rn = 3 are compared to known values 
for ideal states, and that with the closest Euclidean 
distance is taken to be the state. The values for 
pure asters and nematic phases are easy to derive; for 
an aster p(6') = (- cos^, - sin6'), (Qo, Qi, 02, Qa) = 
(0,1,0,0), whereas for the nematic state where p(6') = 
(0,0), (Oo,Qi,Q2,Q3) - (0,0,0,0). Note that while 
an isotropic state would give the same Q^ as for 
the nematic, such states only arise for A:a and P 



well below the considered ranges. For spindles and 
semi-aster states there is a degree of choice in how 
the target Qm are calculated, so we choose sim- 
ple forms that permit exact evaluation of the Qm- 
For the spindle, p(6') — (— cos6',sin0) for 9 G 
(— 7r/4, 7r/4) or 9 E (37r/4, 57r/4) and zero otherwise, 
for which (Qo,Qi,Q2,03) = (0,2/7r2 + 1/2, 0, 2/7r2). 
For the semi-aster, p{9) = {— cos[9/3\, — sm[9/3]) for 
9 e (— 37r/4, 37r/4) and zero otherwise, for which 
{Qo,Qi,Q2,Q3) = (9/7r2,9/47r2,333/17157r2,9/1007r2). 
Variations in these forms have been tested and although 
the boundaries between states shift slightly, the underly- 
ing trends remain the same. 

The occurrence of the four steady-states, plus a fifth 
'vortex' state to be discussed below, with motor density 
and speed are presented in Fig. [3] for two different pres- 
sures. The observed configurations correlate with the 
density and distribution of motors along the filaments. 
The mean number of motors per filament n,^ot/N is plot- 
ted in Fig. [4] and shows an expected increase with the 
attachment rate ^a as well as the pressure. The in- 
crease with pressure can be understood as due to the 
closer packing of the filaments, increasing the number of 
potential attachment points for motors and hence nmot- 
The approximate scaling n-^^t ~ ^a ^^ faster than the 
linear relationship measured for constant volume, two- 
dimensional simulations |19j . presumably due to similar 
reasons: As /ca increases so does the motor density which, 
in this constant-pressure ensemble, allows the system to 
contract, presenting more potential attachment points 
between monomers and hence further increasing the mo- 
tor density. A derivation of the value 3/2 of the exponent 
is not available so far. 

Motors move to the [-|~]-end and dwell there until de- 
taching, thus a greater fraction are expected to occupy 
filament [-t-]-ends, potentially resulting in tight binding 
mediated via many motors. Plotted in Fig. [5] is the frac- 
tion of motors with at least one head at a filament's [-f]- 

end, ninot/'^mot, for the same runs as in Fig. |4J By com- 
paring to the configuration plots in Fig. [3] it is possible to 
infer signatures of crossovers between states in the inflec- 
tion points in these curves. For the slower motors, there 

is an increase in n^Jj/nmot as the nematic state changes 
to a state with a greater degree of polar ordering (spindle 
or semi-aster depending on the pressure). For the faster 

motors there is a marked increase in nJ^Qj/rimot with /ca, 
which corresponds to the crossover to the aster state with 
a high degree of [-|-]-end binding. Further indication of 
the importance of [-|-]-end binding is presented in sec- 
tion IV where enhanced end-unbinding rates /ce > fco 
arc considered. 



B. Dynamics and vortices 

The stationary states described above admit no spon- 
taneous non-equilibrium flows, despite the motor motion 




FIG. 2: Snapshots of steady-states for a low motor attach- 
ment rate kA/kn ~ 1 (left) and a higher rate kx/ku ~ 30 
(right). Conversely, the top line is for fast motors kuTt = 
3.75 X 10~^ and the bottom line for motors 10 times slower. 
Filaments are shaded light (dark) towards their plus (minus) 
ends, respectively. These states are referred to as (a) spindle, 
(b) aster, (c) nematic and (d) semi-aster. The other param- 
eters are P/Pq = 0.03 and ^e = ^d. Movies of the same 
parameter values are available from the supplementary infor- 
mation. 



generating a positive energy flux: The increase in the 
stored motor elastic energy due to motor motion and 
thermal drift of the connected filaments is balanced by 
the loss due to detachment, with no observed net translo- 
cation or rotation of the filaments in steady-state. Col- 
lective rotation of all filaments about a fixed center arises 
for one region of the considered parameter ranges, but ap- 
pears to be a transient flow that irreversibly contracts to 
a non-rotating semi-aster configuration. These states are 
referred to here as vortices due to their superficial simi- 
larity with the rotational modes observed in microtubule- 
kinesin assays |15[ I16j and are described in detail in this 
section. A snapshot is given in Fig. |6] and movies are 
provided as part of the supplementary information [29] . 

Collective rotation of the whole system can be quan- 
tified by the mean angular velocity of filament centre- 
of-mass vectors r relative to the system center, or alter- 
natively by the net transverse velocity of each filament's 
centre-of-mass relative to its polarity, ((v x p)z). Here 
we employ the latter as it is available for all of our runs, 
but we have confirmed that it closely tracks the angu- 
lar velocity in those runs for which both were measured. 
Examples of v x p for 3 independent runs are given in 
Fig. [71 and show finite rotation of either sign until the 
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FIG. 3: Variation of steady-state with motor density and 
speed for (a) P/Po = 0.01 and (b) P/Po = 0.024. Mark- 
ers denote actual states determined as described in the text 
and the boundaries are placed midway between points. The 
vortex region in (b) is a transient configuration that is ex- 
plained in section [III B| and is delineated as those states with 
a vorticity V exceeding 0.7. Since it eventually contracts to a 
semi-aster, the two distinct semi-aster regions in (b) become 
contiguous in steady-state. 



system irreversibly contracts to a semi-aster state and 
rotation ceases. This contraction time can be confirmed 
by visual observation of system states, and can be pre- 
cisely located by fitting the system radius as a function 
of time, R{t), to the four-parameter hyperbolic tangent 
R{t) = i?„i„ + Ai?tanh[(t - t™"*)/Ai]. The mean of 
((v X p)^) is presented in Fig. [Tlfor each run as a hori- 
zontal line segment, that extends from t = to the con- 
traction time f^""' found from this fit. In all cases, i'^""* 
coincides with the rapid decay of v x p to zero, providing 
independent confirmation that rotation ceases when the 
system contracts to a semi-aster. 

It is now possible to define a vorticity order parameter 
for each point in parameter space. For each run a, the 
mean /i^ and standard deviation aa of (v x p)|z is calcu- 
lated starting from t = up to the time that the system 
contracts. This is regarded as significant if the mean is 
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FIG. 4: No. of motors per filament versus attachment rate 
kA/kn for the motor speeds and external pressure denoted in 
the legend. The thick black dashed line has a slope of 3/2. 
Where data for the required P was not available, rimot was 
interpolated from runs with P slightly higher and lower than 
the target value. 



comparable to or larger than the standard deviation, but 
since the sign is arbitrary we also take the absolute value 
to give the vorticity for a single run. 



Va- 



>0. 



(3) 



This is then averaged over all runs with the same pa- 
rameters to give the mean vorticity V = Va- A given 
point in parameter space is then regarded as exhibiting 
a (transient) vortex if V exceeds some arbitrarily-chosen 
value of order unity. The corresponding region of param- 
eter space for V > 0.7 is plotted in Fig. [3] and arises for 
higher densities of motors that arc not so fast that they 
aggregate at [-|-]-ends, which would stabilize an aster rel- 
ative to a vortex. Independently varying the fraction of 
[-t-]-ended motors by increasing the end-detachment rate 
fcE supports the existence of a critical fraction for vortex 



formation, as discussed in Sec. IV 



The reciprocal relationship between vorticity and con- 
traction time is clearly evident when both quantities are 
plotted together; see Fig. [8] Stronger vortices have a 
shorter lifetime than weaker vortices. The distribution 
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FIG. 6: Snapshot of a vortex rotating in the anti-clockwise di- 
rection as presented, for parameters P/ Pq = 0.024, kp^/kY) = 
35 and kyiTb — 7.5 x 0^"^ taken at a time t/rh ~ 5.3 x 10'^. 
The color code is the same as in Fig. [2] 



FIG. 5: Fraction of motors with at least one head at a fila- 
ment's [-|-]-end versus fcA/fco for the same data as in Fig. |4] 
so the top two lines correspond to fast motors and the bottom 
to slow motors. 



of contraction times is presented in Fig. ^a) for a sin- 
gle pressure, and although noisy it appears to be roughly 
consistent with an exponential form suggesting that con- 
traction is a Poisson process. To give some measure of 
the goodness-of-fit, the Anderson-Darling statistic for an 
exponential distribution with an unknown mean has a 
significance level of P « 0.2, whereas that for a normal 
distribution of unknown mean and variance has a sig- 
nificance level of P « 0.025 [301 ISI]- Thus we cannot 
discard the hypothesis that the underlying distribution 
is exponential. Attaining even this noisy data consumed 
considerable computing resources and we were unable to 
repeat this procedure for other parameter values. 

Assuming the true distribution is exponential, this 
would suggest that contraction is triggered by sponta- 
neous fluctuations that occur at a constant rate in time. 
From observation of movies of filament arrangements, a 
likely candidate is the transient void formation frequently 
observed near the outer wall, where nearby filaments are 
attached purely by motors at their [-|-]-ends and not along 
their length. Such voids, when large enough, lead to a 
'hinge'-like mechanism in which the void expands and 
one section of the polarity field inverts, leading to the 
semi- aster. 



The onset of vorticity is also evident in the histogram 
of the signed vorticity, i.e. the Vq. before taking the mod- 
ulus in eqn. (|3|, which can be positive or negative de- 
pending on the direction of rotation. For low pressures 
with V <C 1 this distribution is unimodal around the ori- 
gin, but becomes bimodal when vorticity is more evident 
as demonstrated in Fig. ^Vo) . Of the 20 runs presented 
here, 12 rotated in one direction and 8 in the other, which 
has a significance interval of P « 0.5 as determined from 
a Binomial test with equal probabilities for both direc- 
tions. This is to be expected given our use of stochastic 
initial condition that does not predispose the system to 
any preferred rotational direction. 

Independent confirmation of vorticity can be inferred 
from the mean-squared angular deviation ((A6')^) al- 
ready defined in Sec. |ll] This is plotted in Fig. [lO] for 
the same parameters as above as a function of the lag 
time At, averaged over all waiting times t up until i'^""*. 
There is a crossover from linear behavior (A0)^ ~ At for 
low pressures with low vorticity, to a more rapid scaling 
[AOY ~ (^i)^ for pressures well into the vortex regime. 
Since this quantity is the angular analogue of the mean 
squared displacement for translation degrees of freedom, 
these two limits can be regarded as diffusive and ballis- 
tic, respectively. Microscopically the diffusive limit corre- 
sponds to fluctuations with no net drift, whereas the bal- 
listic limit arises when all filaments are rotating around 
the system with a constant angular velocity in the same 
direction. Therefore the vortex state should correlate 
with ballistic motion, and comparison of Figs. [5] and [TO] 
confirms this. This figure also demonstrates that the in- 
tegrated angular rotation of the vortex before contraction 
is typically larger than 7r/2, much larger than the diffu- 
sive drift w 7r/10 over the same time frame. 
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FIG. 7: Examples of the vorticity (v x p)!^ scaled by a char- 
acteristic velocity v for k^/kiy = 35, kyiTb = 7.5 x 10~^ and 
P/ Pq = 0.024. For clarity error bars are only given for a sin- 
gle run. The thick horizontal line segments denote the mean 
value up to the time when the vortex contracts to a semi-aster, 
and are plotted up to this time. 



FIG. 8: Contraction time to a semi-aster t'^°°'/ri, (left axis; 
solid circles) and vorticity V (right axis; open squares) versus 
pressure P/Po for fcA/^D = 35 and kyiTb = 7.5 x 10"'^. The 
contraction time was determined by the fit of the radius to a 
hyperbolic tangent, or assigned the maximum value of t'^""* = 
1.6 X 10*Ti, if no contraction had occurred within this time. 
Each point represents 5 independent runs. 



IV. ENHANCED DETACHMENT FROM ENDS 



In the simulations that accompanied the microtubule 
experiments, it was claimed that the residence time at the 
microtubule [-f ]-ends played a crucial role in determining 
the vortex stability, with an enhanced end- detachment 
rate required to form vortices [H]. By contrast, for our 
model it is the fraction of motors at filament [-|-]-ends 
that determines vortex stability relative to an aster. The 
vorticity, motor density and fraction of motors at [-1-]- 
ends are plotted in Fig. [TT] against end-detachment rates 
^E > ^D for two sets of /ca, ^m and P. It is clear that 
increasing feg can both destroy a vortex that existed when 
^E = ^D, and create a vortex when k-^ = fco gave an 
aster. In order of increasing /ce, the sequence aster — >■ 
vortex — > semi- aster (where any vortex is either absent 
or too short lived to be discerned) is typically observed, 
although we do not claim this sequence is followed by all 
points in parameter space. There is a slight increase in 
motor density in the semi-aster state as evident from the 
figure, resulting from an increase in potential attachment 
points due to the increased density. 
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FIG. 9: (a) Probability distribution of contraction times on 
log-linear axes for P/ Po = 0.024, kx/ku ~ 35 and kuTt = 
7.5 X 10"^. The thick line gives the best fit to an exponen- 
tial distribution which has a mean « 1.7 x 10*, correspond- 
ing to ~ 0.75 full rotations (the longest vortex survived for 
~ 2.6 rotations). Data corresponds to 20 independent runs, 
(b) Normalised probability histogram of signed vorticity for 
P/Pq = 0.020 (white bars in the background; 5 runs) and 
P/Po = 0.024 (shaded bars in the foreground; 20 runs). 
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FIG. 10: Mean-squared changes in angle {(A6)^) versus lag 
time At for the pressures given in the legend, kp^/k-Q = 35 
and kyifb = 7.5 x 10"'^. The short thick line segments have 
the slope given. 



Thus residence time at [+]-ends, which is oc /cg , is 
not the determining factor with regards vortex stability 
here. Rather, vortices coincide with around 25% of mo- 
tors at [-|-]-ends as highlighted in the figure. There is 
no critical dependency on the motor density, although 
we speculate that below some minimum value spindles 
or nematic states would be observed instead. The crit- 
ical fraction 25% will likely depend on parameters that 
were not varied in this work, such as filament length M 
and the motor spring stiffness. A systematic survey of 
these parameters, or of k^ ^ fcn for all k^, kyi and P, is 
however beyond the scope of this work. 



V. CONTROLLED VOLUME 

One message from the previous sections is that the ob- 
served steady-state is predominately determined by the 
density of motors and the fraction at [-(-]-ends. It may ap- 
pear that the primary role of motor motion, which would 
be the source of any non-equilibrium effects in this model, 
is merely to select the fraction of [-|-]-ended motors, faster 
motors giving a higher fraction. It might even be specu- 
lated that even the transient vortex state is driven, not 



FIG. 11: Variation of vorticity, motor density and fraction of 
[+]-end motors with /ce for /cA/fcD = 35 and kyiTb = 7.5 x 10~^ 
(solid lines, filled squares) and kA/ku = 60 and kuTb ~ 37.5 x 
10"'^ (dashed lines, open diamonds). P/Po = 0.024 in both 
cases. The thick dashed line in the lower plot corresponds to 
25%. Quantities were measured just prior to contraction, or 
in steady-state if there was no contraction or it happened too 
rapidly to discern. 



by motor motion, but rather as a protracted buckling 
event powered by the pressurised walls. 

It is straightforward to show that motor motion can 
drive vortex motion, however. Plotted in Fig. [12] is the 
rotational velocity (vxp)/w for two independent runs in a 
box with fixed radius, where v is a characteristic filament 
velocity. The system is initially in an aster configuration, 
but when the radius is suddenly reduced by 6/2 at a 
time t/Tb « 1.6 X 10^, the system switches to a rotating 
vortex state that appears to be long-lived; the total time 
window in this figure is an order of magnitude longer than 
the longest vortex described in Sec. Ill B (which has the 
same parameters). Since there is no energy input from 
the walls, the only possible cause for this rotation is the 
motor motion. Thus the pressure ensemble is important 
to let the system adjust its density to the vortex state; 
however, the same pressure also destabilizes the vortex 
state, because it favors further contraction into the semi- 
aster. 

Although the magnitude of the rotational velocity re- 
mains fixed (note that the characteristic velocity v is the 
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FIG. 12: (a) Filament rotation for two independent runs at 
fixed volume. The imposed radius R is given in the lower 
panel. Snapshots for the run corresponding to the solid black 
line at points (b) and (c) are given in the lower figures. A 
short time at the initial radius R/b = 38 was required to 
avoid numerical instabilities. 



same for both runs and constant in time), the direction 
aperiodically reverses as evident in the changes of sign 
in the figure. The statistics of time intervals between di- 
rection switching suggest that the underlying mechanism 
may be the same as for the contraction to the semi-aster 
state in the constant pressure case. Specifically, the mean 
switching time Ai™'''=i^/rb = 21.6 x lO-^iS.S x 10"^ is con- 
sistent with the mean contraction time « 1.7 x lO** mea- 
sured earlier, and again is consistent with an exponential 
distribution (significance level P ^ 0.2 from n — 7 values 
using the Anderson-Darling test [3DJ[3T]). This suggests 
that the same spontaneous fluctuation that permits con- 
traction under constant pressure instead promotes rota- 
tional reversal under constant volume. 



VI. DISCUSSION 

It has been demonstrated that the vortices described 
here involve the collective rotation of the filaments about 



a fixed centre. This was predicted by the nematody- 
namics theory of Kruse et al. [17], but contrasts with 
the simulations of intersecting filaments jl6| for which 
movies appear to show no actual filament rotation, rather 
the motors run along a static vortex configuration of the 
growing filaments. It is possible that our inclusion of ex- 
cluded volume interactions, which are necessary to gener- 
ate the nematic elasticity required by the theory but were 
absent in these earlier simulations, may explain this dis- 
crepancy. Those simulations also employed growing fila- 
ments, whereas here as in the continuum modeling these 
lengths were fixed, suggesting a further potential source 
of discrepency. Unlike the continuum theories, however, 
our vortices are only one filament in radius so there is no 
radial gradient in polarity field, making direct compari- 
son problematic. We conclude that the challenging task 
remains to demonstrate a definitive link between macro- 
scopic vortices and microscopic filament-motor interac- 
tions. An important aspect could be the system size of 
the microscopic models, because a minimum size much 
larger than currently accessible by simulations may be 
needed to see vortices in a bulk system. Further exper- 
iments might elucidate the underlying mechanism. For 
instance, fiuorescently tagging a small fraction of micro- 
tubules would allow individual filament rotation (if any) 
to be visualised. 

The variation of steady-state structure with motor 
speed and density shown in Fig. [3] can in large part be 
understood as due to a non-uniform density of motors 
along the filament, in particular the fraction of motors 
at a [-|-]-end which dwell there before detaching. Varying 
the end-detachment rate /cg confirms that a strong bind- 
ing at [-|-]-ends can stabilise an aster relative to a vortex 
or semi-aster. The motor speed /cm plays a role in select- 
ing the distribution of motors along the filament, but also 
contributes to the rotation of vortices as inferred from the 
fixed volume system in Sec. [V] Therefore we claim that 
the observed vortex is a genuine non-equilibrium state 
powered at least partially by the unidirectional motion 
of energy-consuming motor heads along the filaments, 
although at constant pressure they appear to be tran- 
sient. It is not clear if varying some other parameters 
may produce stable vortices. 

Systematically quantifying the role of all of the 
model parameters is clearly challenging for such a high- 
dimensional parameter space, and here we have adopted 
the pragmatic approach of holding most parameters fixed 
while varying those deemed most likely to be critical. 
Eventually the impact of all parameters on structure and 
dynamics will need to be quantified if a broad descrip- 
tion of active gels is to be attained. Here, we highlight 
two parameters likely to reveal novel or interesting be- 
haviour. First, the filament length L — Mb was fixed 
at M — 30 monomers throughout, whereas extensive 
simulations with M — 25 revealed similar steady-state 
diagrams as Fig. p] but no vortices j32j . Increasing the 
aspect ratio therefore seems to enhance vorticity, and it 
would be interesting to quantify this effect. Secondly, 
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the elastic parameters for the wall were set to maintain a 
roughly circular shape, as in the emulsion experiments of 
Pinot et al. |11| . However, in those experiments flexible 
vesicles were also considered that produced a richer ar- 
ray of observed structures, and this effect could be easily 
investigated within our model by lowering the bending 
stiffness of the wall. 
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